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ABSTRACT 


As  part  of  its  continuing  support  of  IMSD 
projects,  the  Aerodynamics  Department  (53-12) 
has  investigated  the  aerodynamic  character¬ 
istics  of  basic  shapes  in  rarified  flows.  This 
report  contains  an  analysis  of  the  drag  and 
moment  forces  on  spheres,  conical  frustums, 
and  cylinders  for  arbitrary  angles  of  attack. 

The  analysis  is  based  on  the  assumption  of 
Newtonian  flow  with  diffuse  particle  reflection. 
Explicit  expressions  are  derived  fc-  quantities 
of  interest. 
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SECTION  1 
INTRODUCTION 


At  altitudes  of  about  75  miles  and  higher,  the  mean  free  path  between 
collisions  of  atmospheric  particles  is  so  large  that  the  concepts  of  con¬ 
tinuum  fluid  mechanics  are  usually  not  valid.  The  computation  of  aero¬ 
dynamic  forces  and  moments  must  then  be  based  upon  concepts  from  the 
kinetic  theory  of  gases,  which  provides  some  understanding  of  the  molecular 
motion  of  the  gas.  The  gas  particles  which  make  up  the  atmosphere  can  be 
considered  in  thermal  equilibrium  in  a  given  region,  and  the  individual 
particles  possess  velocities  which  are  described  by  the  Maxwell-Boltzmann 
velocity  distribution  function.  This  function  makes  it  possible  to  com¬ 
pute  the  number  of  particles  that  strike  a  surface  of  unit  area  per  unit 
of  time  and  the  resulting  momentum  transfer  if  the  nature  of  the  molecular 
reflection  process  at  the  surface  is  known.  It  is  generally  assumed  that 
particles  reflect  from  surfaces  either  specularly  or  diffusely,  the  former 
meaning  that  a  particle  has  an  equal  angle  of  incidence  and  angle  of  re¬ 
flection,  while  the  latter  implies  that  the  reflected  particles  leave  the 
surface  randomly.  A  particle  undergoing  specular  reflection  produces  a 
momentum  exchange  which  is  proportional  to  twice  the  normal  component  of 
velocity,  and  the  computations  are  relatively  simple  For  diffuse  reflec¬ 
tion,  it  is  assumed  that  an  incident  particle  hits  the  surface  and  comes  to 
rest,  giving  up  its  total  momentum,  and  then  is  ejected  in  a  manner  that  is 
related  to  the  surface  temperature.  Diffuse  reflection  is,  of  course,  a 
more  complicated  process  than  specular  reflection  and,  unfortunately,  ex¬ 
perimental  results  pertaining  to  flight  conditions  m  the  upper  atmosphere 
are  not  available.  In  reality,  both  types  of  reflection  may  occur,  although 
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The  diffuse  is  probably  predominant.  Until  more  data  are  available,  it 
seems  reasonable  to  assume  that  the  reflection  process  is  completely 
diffuse.  This  means  that  the  particles  which  are  ejected  from  a  surface 
possess  a  Maxwellian  distribution  corresponding  to  the  surface  temperature. 

Free  molecular  flow  consists  of  superimposing  a  uniform  flow  on  the  thermal 
motion  of  the  particles  and  using  the  superimposed  distribution  function  to 
calculate  the  forces  arising  from  the  momentum  exchange  due  to  collision  of 
particles  with  the  surface  of  tbe  vehicle,  Boundary  layers  and  shock  waves 
do  not  exist  in  a  clearly  distinguishable  form,  because  the  mean  free  path 
is  quite  long  at  high  altitude;  hence,  all  intermolecular  collisions  are 
neglected. 


In  the  kinetic  theory  of  gases,  a  quantity  called  the  most  probable  molecular 
velocity  is  defined.  This  quantity  (C  =  J  2RT),  a  sort  of  average  velocity, 
is  proportional  to  the  square  root  of  the  temperature,  and  is  related  to  the 
speed  of  sound.  The  molecular  speed  ratio  is  defined  as  the  uniform  flow 
velocity  divided  by  the  most  probable  molecular  velocity  and  is,  of  course, 
related  to  the  Mach  number .  For  particular  geometries  such  as  the  cone, 
cylinder,  and  sphere,  the  drag  coefficient  as  determined  from  free  molecular 
theory  is  a  function  only  of  the  molecular  speed  ratio,  assuming  that  the 
vehicle  surface  has  a  given  temperature.  For  example,  the  drag  coefficient 
for  a  sphere  in  free  molecular  flow  with  diffuse  reflection  is  given  in 
Reference  1*  as  follows’ 


C 


D 


+  Fi? 


+  2  (i  +  s?  -  riir  ) 


erf  S  + 


2 jr 

5Sr 


(1) 


In  Equation  (l),  S  is  the  molecular  speed  ratio,  and  is  a  molecular  speed 
ratio  based  on  the  wall  temperature  and  refers  to  the  reflected  particles  only. 
Assuming  that  only  altitudes  of  over  100  miles  are  of  interest,  it  is  known 

*  See  References 


1-2 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


LMSP-513^ 


(Reference  2)  that  the  surface  temperature  of  a  vehicle  is  determined 
almost  entirely  by  radiation  phenomena.  From  the  definition  of  S  and  S^, 
the  following  is  obtained: 


where  is  the  surface  temperature,  and  1\  is  the  ambient  temperature. 
Assuming  that  the  surface  emissivity  is  such  that  T^  =  >yC  R  and  noting 
from  the  ARDC  atmospheric  tables  that  T  is  12D0°R  or  larger  above  100 
miles,  Equation  (2),  using  T  =  1200°R,  becomes 

S  =  1.^7  S  (3) 

r 

and  Equation  (l)  can  be  plotted  as  shown  in  figure  1  by  the  upper  curve. 
This  curve  is  a  maximum  for  altitudes  above  100  miles,  since  a  higher  T^ 
would  cause  S  to  increase,  thus  reducing  the  contribution  of  the  last  tein 
in  Equation  (l)  to  the  drag  coefficient.  Taking  the  limit  of  Equation  (l) 
as  S  approaches  infinity  results  in  =  2.  It  is  apparent  that  the  last 
term  decreases  more  slowly  than  the  others  because  it  is  proportional  to  S 
while  the  others  are  exponential  or  of  higher  exponent  The  last  term  in 
Equation  (l)  arises  from  consideration  of  the  diffuse  reflection,  and  the 
above  simplifications  applied  to  that  term  are  valid  only  when  the  surface 
temperature  is  primarily  controlled  by  radiation. 

Although  a  free  molecular  analysis  can  be  applied  to  any  simple  shape,  some 
times  the  results  cannot  be  written  in  terns  of  elementary  functions,  and 
numerical  methods  must  be  used.  This  situation  occurs  in  particular  for  a 
cone  at  an  arbitrary  angle  of  attack  It  is  c  sirablo  then  to  have  an 
approximate  method  of  determining  forces  and  moments  on  vehicles  of  various 
shapes.  It  is  well  known  that  the  Newtonian  flow  concept,  assuming  that 
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0  2  4  6  5  10  12  14  16 
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Figure  1  Drag  Coefficient  vs  Molecular  Speed  Ratio 
for  a  Sphere 

the  particles  come  to  rest  on  the  surface,  produces  a  drag  coefficient  of 
2  for  all  bodies  based  on  the  projected  frontal  area.  This  corresponds  to 
the  limiting  case  of  Equation  (l)  as  expected,  because  thermal  motion  of 
the  particles  is  neglected,  and  hence  can  serve  as  an  approximation  as  long 
as  S  is  relatively  large.  For  example,  when  S  =  10,  using  CD  =  2  provides 
a  result  which  is  4  percent  low,  as  seen  from  tne  upper  curve  of  figure  1. 

So  far  nothing  has  been  stated  concerning  the  reflection  of  particles  in  New¬ 
tonian  flow.  In  the  Newtonian  hypersonic  approximation,  it  is  assumed  that  the 
particles  give  up  their  normal  component  of  velocity  and  then  slide  off  the  sur¬ 
face  tangentially.  Newtonian  flow  with  specular  reflection  imparts  twice  as  muc 
momentum  to  the  surface  as  the  Newtonian  hypersonic  approximation.  Newtonian 
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flow  with  completely  diffuse  reflection  is  more  realistic  for  the  problem 
under  consideration,  and  the  sphere  will  be  analyzed  in  detail,  below  us  an 
illustration  of  the  method. 

The  total  momentum  transferred  to  an  annular  element  is  the  drag  on  that 
element,  and  integrating  over  the  spherical  surface  gives 

V  2 

C  =  ^ — jj - x  =  I  4  cos  0  sin0d0-  2 

U  |pr*r  Jo 


(4) 


where  P  is  the  atmospheric  density,  U  the  velocity,  and  r  the  radius  of  the 
sphere.  For  diffuse  reflection,  the  kinetic  theory  of  gases  provides  an 
expression  for  the  pressure  of  the  ejected  particles,  assuming  that  they  are 
in  thermal  equilibrium  with  the  wall. 


Pr 


m 


(5) 


ihe 

iur- 

much 

i 


where  mr  is  the  mass  leaving  a  unit  surface  area  per  unit  of  time,  assume] 
equal  to  the  incident  mass,  and  =  \J2  RTr  .  Therefore, 


P 


r 


/FPU2 

2  Sr 


cos  0 
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and  taking  the  component  of  thia  term  in  the  drag  direction  and  integrating 
over  the  surface  givea 


D 

r 


PU2  f  /2  2  2 

— g-  I  2  nr  cos 

r  Jo 


6  ain  0d© 


(6) 


or 


The  total  drag  coefficient  ia  then  the  aum  of  Equationa  (4)  and  (6) 


°D 


2  + 


2j£ 

3  Sr 


(7) 


Thia  equation  ia  plotted  in  figure  1,  and  the  curve  ia  labeled  Newtonian- 
diffuae.  The  curve  ia  a  fairly  good  approximation  of  the  free  molecular 
result,  and  even  at  S  =  6  there  is  only  a  2. 8-percent  discrepancy.  Thia  ia 
becauae  the  reflection  terra  doea  not  decreaae  rapidly  aa  S  becomea  large. 

The  utility  of  the  Nevtonian-diffuae  analyaia  can  be  appreciated  by  atudying 
figure  2,  which  ia  a  plot  of  altitude  againat  the  molecular  apeed  ratio 
correaponding  to  the  velocity  of  a  aatellite  in  circular  orbit  at  each  alti¬ 
tude.  Aa  can  be  aeen,  the  Nevtonian-diffuae  approximation  providea  a  result 
for  the  drag  coefficient  of  a  aphere  which  ia  at  moat  about  3  percent  in 
error  throughout  the  altitude  range  in  which  aerodynamic  forcea  on  aatellitea 
are  of  aignificance .  At  the  lower  altitudea,  this  error  becomes  aomew’nat 
leaa,  and  the  Newtonisn-diffuse  theory  ahould  be  adequate  for  moat  purpoaea 
even  when  the  aerodynamic  moment  ia  a  critical  factor.  The  error  aaaociated 
with  the  Nevtoniaa-Diffuae  theory  ia  different  for  other  bodiea  and  ia  dia- 
cuaaed  ia  the  Appendix. 
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Figure  2  Molecular  Speed  Ratio  vs  Altitude 


The  following  sections  provide  formulas  for  computing  the  aerodynamic 
moment  and  the  drag  forces  on  shapes  which  usually  are  of  interest  for 
orbiting  vehicles. 
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SECTION  II 

ANALYSES  FOR  CONE  SHAPE 


The  analysis  of  drag  and  moment  on  a  cone  is  divided  into  two  parts,  one 
for  the  angle  of  attack  greater  than  the  semi -vertex  angle,  and  one  for  the 
angle  of  attack  less  than  the  semi-vertex  angle.  The  moments  due  to  the 
incidence  and  reflection  of  particles  from  the  surface  are  also  considered 
separately,  and  the  results  given  below. 

2.1  INCIDENT  MOMENT  Ct  <  5 

Figure  3  illustrates  the  geometry  of  the  problem  and  the  notation  used  in 
the  following  equations.  The  force  on  a  unit  of  surface  area  is  in  the  free 
stream  direction  and  is  given  by 

"l"  j  x  (6) 

where  i  is  the  direction  cosine  be twee i  the  free  stream  and  the  inward  normal 
x 

at  any  point  on  the  cone: 

jP  =  sin  i  cos  a  +  si  n  ct  cos  a  cos  *  (9) 

X 

The  moment  arm  from  this  force  to  an  arbitrary  vehicle  center  of  gravity 
located  on  the  cone  axis  a  distance  L  from  the  vertex  is 

(L  -  y)  sin  a  -  y  tan  6  cos  a  cos  <t>  (10) 
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Figure  3  Geometry  and  Notation  for  Cone 


and  an  element  of  surface  area  is 


ydy  d* 


n) 


Combining  the  above  four  quantities  gives 


mcg  =  2p 


^  t<L " y)  sln  “ ' y  tan  *  oos  °  cos*  J 

0  yl 


tan  i 
cos  a 


y  dy  d$ 


(12) 
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and  carrying  out  the  indicated  integration  provides  the  following  result: 


M^g  -  2pU^  it  t&nd  6  sin  a  cos  a 


m 


(13) 


where  y^  £  0  for  a  cone  frustum. 

2.2  DIFFUSE  REFLECTION  M0MEI7T  a<  6 


The  reflection  process  exerts  a  pressure  on  the  surface  which  it;  given  by 


where  the  incident  mass  is  equal  to 


m 


i 


-  P  \j£ 

x 


(14) 


05) 


The  pressure  is  divided  into  components  parallel  and  norma],  to  the  cone 
axis,  giving  a  combined  moment  arm  equal  to 

(16) 


r  .  ^  y  COS0 

L  cos  6  cos  0  -  * - 


and  the  moment  Is  expressed  as  follows: 


-*_) 

cos  6  / 


cos  <t> 


tan  6 
cos  6 


y  dy  (1  <P 


07) 
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The  result  of  the  integration  is 


n  tan  < 
2  cos  s 


sin  a 


T  2  , 

L  cos  6 


2-3  INCIDENT  MOMENT  <*>$ 

The  results  can  he  obtained  for  this  case  in  a  manner  similar  to  that  for  the 
case  a  <  i  except  that  now  part  of  the  cone  is  shielded  from  the  free  stream. 
This  area  is  defined  by  the  locus  of  points  at  which  the  free  stream  becomes 
tangent  to  the  cone  and  is  found  by  setting  Equation  (9)  equal  to  zero,  solving 
for  and  using  this  as  the  limit  of  the  $  integration.  The  final  result  for 

this  case  is 


tan  6 
cos  6 


[ 


A  6in  a  +  B  cos  a  tan  t 


(19) 


where  A  =  cos  a  sin  i 


IT  -  C08~^ 


) 


(cot  a  tan  i  )  +  sin  a  cos  6  sin  cos"  (cot  atari  6  ) 


c 


) 


and  B  =  cos  a  sin  i  sin 

+  sin  a  cos  j 


£  cos  1  (cot  a  tan  i  )J 

^  i  cos"1  (cot  a  tan  i  )  -  ^  sin  £2  cos-1  (cotatanfljj 
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2.4  DIFFUSE  REFLECTION  MOMENT  Ct 

The  moment  about  an  arbitrary  point  located  on  the  line  of  the  cone  axis  is 


where  B  is  defined  above. 

2.5  DRAG  FORCE  d<0 

The  drag  force  is  made  up  of  contributions  from  the  incident  and  reflected 
particles.  The  result  below  includes  both  contributions. 


If  the  drag  coefficient  is  desired,  Equation  (2l)  must  be  divided  by  —  PU^  A, 
where  A  is  some  particular  reference  area.  Usually  it  is  most  convenient  to 
use  the  base  area  for  conical  shapes,  and  hence  for  this  case 

A  =  n(y2  "  yl)  ta”2  4 
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2.6  DRAG  FORCE  a  >  5 


The  drag  force  for  conditions  when  a  >  6,  including  the  incident  and  diffuse 
contributions,  is  given  below. 


D  =  2pU 


^2  "  y2l^ 


tan  S 
cos  5 


cos  a  sin  & 


^  it  -  cos  ^  (cot  a  tan  5)J 
+  I  —  cos  a 


VTi 


[ 


-1 


+  sin  a  cos  6  sin  cos  (cot  a  tan  6) 


1+2  cos  a  sin  &  — 

CO 


sin  5 


5T  (22) 


4.  fji  —2  ___2e  |  It  l  ___-l,  .  .V  1  .  L  -l 


2S/^T  sin^  a  cos"~&  j  §  "  f  cos  X(cot  01  tan  5)  -  j-  sin  |2  cos  '"'(cot  a:  tan  &  )] 
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SECTION  III 

ANALYSES  FOR  CYLINDRICAL  SHAPE 

The  geometry  and  notation  for  the  cylinder  are  shown  in  figure  k,  and  the 
essential  steps  in  the  derivation  of  the  equations  are  given  below. 

3.1  INCIDENT  MOMENT 

The  direction  cosine  between  the  free  stream  and  an  inward  normal  of  the 
surface  is 


i  =  sin  a  cos  ®  (23) 

X 

The  moment  arm  from  an  arbitrary  point  on  the  cylinder  axis  to  an  element 
of  the  surface  is 


(L  +  y)  sin  a  -  r  cos  a  cos  $ 


(24) 


and  hence  the  moment  can  be  written  as 

y2  ,  */2 

.2  ' 


mcg  ■  2?u 


r  e  '  r 

/  / 

J  Q  U 


y)  sin  a  -  r  cos  a  cos  $ 


sin  a  cos  $  r  d$  dy  (25) 


The  result  of  the  integration  provides  the  following  result: 


MCG  =  2pU  ry2  8ln  a 


7t  ^2^ 

(L  +  —  I  sin  a  -  cos  a 


(26) 


If  it  is  desired  to  compute  the  moment  about  some  point  located  inside  the 
cylinder  a  distance  L  from  the  origin  of  the  y  coordinate,  then  a  negative 
sign  is  inserted  in  front  of  L  in  Equation  (26).  The  ends  of  the  cylinder  are 
not  included  in  the  analysis. 


3-1 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  or.d  SPACE  DIVISION 


LM3D-‘j]  32 


Figure  k  Geometry  and  Notation  for  Cylinder 


3.2  DIFFUSE  REFLECTION  MOMENT 


Using  the  definition  of  the  pressure  due  to  reflection  as  given  in  Equation  (l4), 
the  following  equation  results  for  the  moment  about  an  arbitrary  point  on  the 


cylinder  axis: 


M 


CG 


■f 


F7  el  /V 

VTi  s  /  / 

Jo  Jo 


n}2 


(L  +  y)  sin  e»  cos^  #  r  d$dy 


(27) 


or 


ry^  sin  a 

C. 


(20) 


As  before,  the  moment  about  a  point  inside  the  cylinder  is  obtained  by  putt i  tq 
a  negative  sign  in  front  of  L.  No  contribution  from  the  end6  of  the  cy.l  indcr 
is  included  in  the  above  computations. 
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3.3  DRAG  FORCE 


The  drag  force  equation  given  below  consists  of  contributions  from  he 
incident  and  reflected  particles. 


D  = 


sin 


*] 


(29) 


The  ends  of  the  cylinder  are  neglected  in  Equation  (29)>  as  can  be  observed 
since  the  drag  is  zero  at  zero  angle  of  attack. 


3-3 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  ond  SPACE  DIVISION 


LM3D-‘>l.i2 


SECTION  IV 

ANALYSES  FOR  SPHERICAL  SHAPE 


Since  the  net  momentum  exchange  on  a  sphere  produces  only  a  drag  f<  rce, 
those  computations  in  which  an  entire  sphere  is  involved  are  relatively 
simple.  The  drag  coefficient  js  given  by  Equation  (7),  based  on  the  projected 
frontal  area,  and  hence  the  moment  of  the  sphere  about  any  other  point  amounts 
to  determining  the  moment  arm  and  multiplying  by  the  drag. 

In  a  number  of  cases  where  a  vehicle  is  composed  of  a  combination  of  cones, 
cylinders,  or  spheres,  a  cone  or  a  cylinder  may  be  capped  by  a  part  of  a 
sphere  or  a  hemisphere.  In  the  paragraphs  below,  the  method  of  computation 
for  such  cases  will  be  presented. 

4.1  INCIDENT  MOMENT  «<7T/2  - 

Figure  5  illustrates  the  notation  which  is  used.  When  a  frustum  of  a  cone 
is  capped  by  a  portion  of  a  sphere,  the  Bemi-vertex  angle  of  the  cone  defines 
the  angle  8^,  which  determines  the  required  part  of  the  spherical  cap  needed 
to  cover  the  cone.  When  the  angle  of  attack  is  less  than  the  angle  «/2  -  8-^, 
then  the  free  stream  strikes  the  entire  surface  of  the  spherical  cap,  and  the 
following  analysis  is  applicable.  The  direction  cosine  between  the  free  stream 
and  an  inner  normal  is 

£  -  cos  &  cos  a  +  sin  8  sin  a  cos  *  (30) 

X 

and  the  moment  arm  is 

R  =  P  sin  a  +  (r  cos  8  tan  a  -  r  sin  8  cos  $  )  cos  a  (3l) 
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Figure  5  Geometry  and  Notation  for  Sphere 


An  element  of  area  on  the  spherical  surface  can  be  written  a6> 

dA  =  r^  sin<jd0d$  (32) 

and  hence  the  moment  about  an  arbitrary  point  on  the  polar  axis  is 


Carrying  out  the  indicated  integration  provides  the  following  result: 

Mcg  =  2PU^ rrr2  sin  a  cos  a  sin2©^  f  ^  |  cos  (3*0 


4-2 


LOCKHEED  AIRCRAFT  CORPORATION 


MISSILES  and  SPACE  DIVISION 


LMSD-^132 


4.2  DIFFUSE  REFLECTION  MOMENT  Of <77/ 2  -  Q 


Using  the  definition  of  the  pressure.  Equation  (l4),  and  dividing  the  pres¬ 
sure  into  components  perpendicular  and  parallel  to  the  polar  axis  yields  the 
following  result: 


(35) 


Carrying  out  the  integration  gives 


sina 


cos  e 
3 


(36) 


4.3  INCIDENT  MOMENT  Ot  >  7T/2  - 


For  this  case,  the  integration  must  be  done  in  two  parts,  because  over 
part  of  the  body  the  flow  hits  elements  of  area  for  ',11  values  of  $ ,  while 
on  the  rest  of  the  body  the  shielding  effect  does  not  allow  impingement  of 
the  free  stream.  The  moment  is  thu6  given  by 


i  R  dA 

X 


/7T-C0S  ^  (cot  a  coto) 

l  R  dA 
X 


(37) 
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The  integration  is  very  lengthy,  although  it  can  bo  done,  with  the  result 
below. 


sin  a 


^CG  ~  r2  j^p  sin  a  cos  a  +  2 sin  a  cos2a  J 

+  ip  sin  o  cos  a  +  2  sin2  ©J  +  sin^  a 

jin  a  -1  /  >  T  o  ^  o 

“ gp~  cos  a  cos  (cot  a  cot  fl^)  -3*  sin  ©1  +  2r  cos'3  8 

-  3r  cos  91  ^1  -  i  cos2 

-  ^  sin  a  cos2  a  sin-1  £l+2  cot2  a  -  2  esc2  a  sin2 ©^ 

+  ^  sin  a  |  2  sin2a  -  cos41  a  )  ^  (sin2©!  esc2  a  -  cotter) ^ 


-  3r  cos  ©, 


^  sin  a  cos2 


1+2  cot2  a  -  2  esc2  a  sin2fl^J 


2  „ 

_  cos 

sin^a  tan  \ - - - — -  -  j  sin  a  cos 

\lnin  ©^ -cos  a  ’ 


asih-1 


2  2 
2  cos  a  esc  © 


?  /  1-cot2  a  cot2  0,  p 

sin  a  cos 5 - 5 -  +  sin  a  cos  a  sin  0^  cos  0  (38) 

t  cot  a  csc^  0X  d  [  1  1 


-  sin  a  cos  a 
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This  equation  can  be  checked  at  two  points.  When  a  =  n/2  -  0^,  this 
equation  must  reduce  to  Equation  (3^-)  with  oc  =  jt/2  -  0^  in  that  equation 
also.  When  a  =  0^^  =  n/2,  Equation  (38)  must  reduce  to  the  moment  of  a 
hemisphere  about  a  point  a  distance  i  from  the  center  of  the  sphere  which 
is  given  below: 


These  conditions  are  satisfied  by  Equation  (38)  if  the  terms  involving  the 
arc  tangent  are  handled  properly.  In  each  of  those  terms,  the  negative 
root  of  the  radical  is  taken,  and  the  angle  determined  is  a  negative  angle. 
For  example,  if  the  quantity  inside  the  radical  is  one,  then  the  negative 
root  is  minus  one,  and  measuring  the  angle  in  the  negative  direction  gives 
tan  1  (-1)  =  -  n/4.  Also  principal  values  of  the  arc  sine  are  to  be  taken. 

4.4  DIFFUSE  REFLECTION  MOMENT  a  >  n/2  -0 
The  result  for  this  case  is 


Here  again,  in  the  term  involving  arc  tangent,  the  negative  root  is  taken,  and 
the  result  is  taken  as  a  negative  angle. 
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4.5  DRAG  FORCE  a  <  IT/2  -  ^ 

The  drag  force,  including  both  the.  incidence  and  diffuse  reflection  of 
particles  from  the  surface,  is 


4.6  DRAG  FORCE  C* >  7r/2  -  d1 

The  drag  force  for  this  case  is  rather  lengthy,  and  both  the  incident  and 
diffuse  contributions  are  included. 
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D  =  2pU2r2  ^  cos  a  sin2  8.  -  sin2  0(  c<  s~'  (cot  or  cot  0, ) 


n  I  2  .  .  2  \ 

+  g  ^cos  a  +  2  sin  or  f 


ci'G^o  .  -1  r,  „  .2 

— t: -  sin  1+2  cot  a  -2  esc 


‘  asliT  0!  J 


eina 


-  cos  0^  sin2  0^  cec2  a-  cot2< 


sin2  -  .  -1  ,  ^  *1 

+  -an  y—TZ - T 

»  sin  01  -  cos  a 


+  puV  .  f  cos2«  +  f  sin2  a  (2  +  5cos3  8X  ) 


2 _ 3  o  cos"  a  3^-1/ 

^  cos  0^  +  —x -  cos  0^  cos  (cot  a  cot  Q  ) 


i  2 

sin  a 

+  — 2 -  008 


0-^1  -  cos  0^  j  cos”1  (cot  a  cot  0^) 


cos2 a  ,  -1  [  „  2  2  „ 

— 7 sin  -2  cos  a  esc  0 


i-] 


2  , 

cos  a  [ 

- ^ 


2  cos2 a  -  sin2a)\Jkj^£^2iS 

V  cot  or  esc  0, 


1  3  \  I  2  2  2 

+  —  sin'  a  cos  a  V  (sin  61  esc  a  -  cot  a 


)3-! 


sin  '  a  cos  0, 


1  .2  -ill-  cot2  a  cot2  0,. 

+  -  sin  a  tan  \  / _ l_ 

*  cot**  a  esc2  0^ 
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SECTION  V 
CONCLUSIONS 


In  some  of  the  results  for  the  drag  force,  the  contribution  of  the  reflec¬ 
ted  particles  results  in  cumbersome  equations  which  are  difficult  to  evaluate. 
Since  the  drag  is  mainly  of  interest  in  determining  orbit  lifetimes,  it  may  be 
entirely  satisfactory  to  use  the  value  of  2  for  the  drag  coefficient  based  on 
the  projected  frontal  area.  However,  the  effect  of  diffuse  reflection  on  aero¬ 
dynamic  moments  is  not  always  small  in  relation  to  the  total  moment,  and  it 
should  not  be  neglected  if  attitude  control  of  the  vehicle  is  a  critical  de¬ 
sign  factor. 

All  of  the  equations  above  are  valid  for  any  angle  of  attack.  When  they  are 
applied  to  a  particular  composite  body,  however,  the  orientation  of  the  co¬ 
ordinate  axes  must  conform  to  the  system  used  in  the  derivation  of  the  equa¬ 
tions. 
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APPENDIX 


The  equations  presented,  in  this  report  based  on  the  Newtonian-diffuse  theory 
are  of  course  approximate  in  that  the  thermal  motion  of  the  molecules  in  the 
free  stream  is  neglected  This  approximation  introduces  an  error  into  the 
results  which  depends  on  the  tody  shape  and  varies  with  the  angle  of  attack 
for  a  particular  shape.  An  illustration  of  this  can  be  seen  by  inspecting 
the  free  molecular  equations  for  lift  and  drag  on  a  flat  plate  and  observing 
that  the  important  parameter  is  the  product  of  the  molecular  speed  ratio  and 
the  sine  of  the  angle  of  attack.  The  Newtonian-diffuse  equations  for  lift 
and  drag  on  a  flat  plate  can  be  easily  obtained,  and  a  comparison  v.ilh  ‘-he 
free  molecular  equations  indicates  that  the  error  increases  with  decreasing 
angle  of  attack  for  a  given  value  of  S.  Hence, the  resulting  equations  pre¬ 
sented  in  this  report  for  the  cylinder  can  be  expected  to  show  a  similar 
tendency.  Accurate  results  for  a  cylinder  can  be  obtained  by  considering 
free  molecular  theory  which  accounts  for  the  free  stream  thermal  motion  of 
atmospheric  particles.  The  analysis  of  the  aerodynamic  moment  is  presented 
in  the  following  paragraphs,  and  the  result  can  be  expressed  in  terms  of 
Bessel  functions  which  are  easily  evaluated  Applying  a  similar  analysis  to 
the  cone  or  spherical  cap  results  in  expressions  requiring  numerical  integra¬ 
tion.  Therefore, the  Newtonian-diffuse  eiuations  for  those  two  bodies  repre¬ 
sent  the  only  analytical  result  available  for  the  flow  regime  under 
consideration. 

The  notation  and  coordinate  system  used  below  corresponds  to  that  used  by 
Talbot  on  page  U58  in  the  June  1957  issue  of  the  Journal  of  the  Aeronautical 
Sciences.  Talbot's  y  axis  is  used  for  a  local  coordinate  system  and  does 
not  appear  in  the  results;  thus,  no  confusion  will  result  by  using  the  y 
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axis  aa  the  cylinder  axis,as  illustrated  in  Figure  k  of  this  report.  Talbot 
defines  the  N  and  T  directions  to  be  normal  and  along  the  cylinder  axis 
respectively  and  lying  in  the  plane  of  the  pitching  moment,  which  is  the 
quantity  to  be  determined.  Taking  the  T  direction  to  be  opposite  from 
that  of  Talbot,  the  components  of  the  pressure  and  shear  forces  in  the  N 
and  T  directions  can  be  written  in  terms  of  the  angle  of  attack  P  and  the 
polar  angle  $  as  follows: 


PN  /sin  p  sin  4>  1  /  r\  . 

—  =  ( - *7= -  +  —  '  TJ  6ln  ♦  6 

ipu  V  SVrT  2S  V 


2  2 
sin 


+  sin  <j>  +  sin2  p  sin2  ♦  +  ~  sin  p  sin 

V2S2  2S 


/)[** 


(1A) 

erf  (t  sin  t)l 


TN  sin  p  cos2  <J>  f  -|2  sin2  <j>  r  .  .  .  I\  ^  .  X\1  1 

_ — _  e  - -j  e  +  vn  5  sin  $  1  +  erf  (  fc  sin  <J> )  b 

^pu  SVjT  "  L  J  ^ 


(2A) 


tT  _  cosjf-S2  aln2  *  +  ^ 


\po2  S  -G 


|  sin  1  +  erf  (i  sin  $) 


]} 


(3A) 


where  |  =  S  sin  p  and  p^  =  0.  Hence,  the  pitching  moment  about  a  point  on  the 
cylinder  axis  a  distance  L  from  the  origin  of  the  y  coordinate  is  given  by 
the  following  expression: 


it/2  r 
”cc / 

-n/2  0 


(PN  -  TN)(y  +  L)  -  TTr  sin  0 


r  dy  d$  (  4a) 
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The  integrals  in  (4A)can  be  easily  evaluated  giving  the  final  result: 


=  2iya  vf  +  L )  sin  f> 

-pU  \  /  I  i. 


(5A) 


If 

2  2 
+  3  sin  p  e 


[(f»  +  t)  lo  (t)  +  (3  +  Se)ii(V)]}-"%  8ln  p 


COB 


/2v 

where  I  (-&-  ) 
o  \2  ) 


and  I 


are  modified  Bessel  functions  of  the  first  kind 


and  zero  and  first  order  respectively.  When  |  becomes  very  large  or  when 
terms  of  the  order  if  i  can  be  neglected,  then  (pA)  reduces  to  the  Newtonian- 
diffuse  result  given  by  the  sum  of  equations  (26)  and  (28).  A  comparison  of 
these  two  results  shows  that  by  accounting  for  the  free  stream  thermal  motion, 
the  first  term  of  (26)  goes  into  the  more  complicated  second  terra  of  (pA), 
which  is  a  function  of  S  sin  p,  Therefore, it  is  evident  that  the  Newtonian- 
diffuse  result  will  be  in  error  to  an  increasing  extent  as  the  angle  of 
attack  decreases  at  a  particular  molecular  speed  ratio. 
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